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Abstract 

We present systematic results, based on density functional calculations, for the structure and 
energetics of ^He and ^He nanodroplets doped with alkaline earth atoms. We predict that alkaline 
earth atoms from Mg to Ba go to the center of ^He drops, whereas Ca, Sr, and Ba reside in a 
deep dimple at the surface of ^He drops, and Mg is at their center. For Ca and Sr, the structure 
of the dimples is shown to be very sensitive to the He-alkaline earth pair potentials used in the 
calculations. The 5s5p <— 5s^ transition of strontium atoms attached to helium nanodroplets 
of either isotope has been probed in absorption experiments. The spectra show that strontium 
is solvated inside ^He nanodroplets, supporting the calculations. In the light of our findings, 
we emphasize the relevance of the heavier alkaline earth atoms for analyzing mixed ^He-^He 
nanodroplets, and in particular, we suggest their use to experimentally probe the ^He-^He interface. 
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I. INTRODUCTION 



Optical investigations of impurities in liquid helium have drawn considerable attention 
in the past.- In recent years, experiments involving helium nanodroplets have added new 
input into the interaction of atomic impurities with a superfluid helium environment .-i^ In 
particular, the shifts of the electronic transition lines represent a very useful observable to 
determine the location of the foreign atom attached to a helium drop. 

While most impurities are found to reside in the interior of helium droplets,-'^'^'^ it is 
well-established that alkali atoms, due to their weak interaction with helium, reside in a 
'dimple' at the surface of the drop for both helium isotopes.-i2i>i^ The question of solvation 
versus surface location for an impurity atom in liquid He can be addressed within the 
model of Ref. [ill, where a simple criterion has been proposed to decide whether surface or 
solvated states are energetically favored. An adimensional parameter A can be defined in 
terms of the impurity- He potential well depth e and the minimum position Tmin, namely, 
A = permm/(2^^^o"), where p and a are the density and surface tension of bulk hquid He, 
respectively. The threshold for solvation in 4He A ~ Ao, with Ao = 1.9. When A < Aq, 
a stable state of the impurity on the droplet surface is expected, whereas when A > Aq, the 
impurity is likely to be solvated in the interior of the droplet. Impurities such as neutral alkali 
atoms, that weakly interact with helium, are characterized by values of A much smaller than 
the above threshold; their stable state is thus expected to be on the surface of the droplet, 
as experimentally found. 

The shape of the impurity- He interaction potential, however, is not given consideration by 
this model. For cases in which the value of A does not lie near (say, within 0.5) the solvation 
threshold Aq, the shape of the potential surface does not need to be taken into account, 
as the model is predictive outside of this threshold window. However, for values which lie 
close to Ao, consideration of the shape of the potential energy surface, as well as the well 
depth and equilibrium internuclear distance, is mandatory, and more detailed calculations 
are needed to ascertain whether the impurity is solvated or not. It is worth noticing that 
the above criterion works for either helium isotope, although so far, it has been applied to 
^He because experimental data for ^He only appeared recently.-'i^'i^i^ 

Among simple atomic impurities, alkaline earth (Ake) atoms play a unique role. While, 
for example, all alkali atoms reside on the surface and all noble gas atoms reside in the 
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interior of drops made of either isotopej-i^ the absorption spectra of heavy alkahne earth 
atoms Ca, Sr, and Ba attached to a ^He cluster clearly support an outside location of Ca 
and Sr~ and likely also of Ba/^ whereas for the lighter Mg atom, the experimental evidence 
shows that it resides in the interior of the ^He droplets.— 

According to the magnitude of the observed shifts, the dimple in the case of alkaline earth 
atoms is thought to be more pronounced than in the case of alkalis, indicating that alkaline 
earth atoms reside deeper inside the drop than alkali atoms. This will be corroborated by 
density functional calculations presented in the Theoretical Results. Laser-induced fluores- 
cence results for Ca atoms in liquid ^He and ^He have been recently reportecP^ and have 
been analyzed using a vibrating 'bubble model' and fairly old Ca-He pair potentials based 
on pseudopotential SCF/ CI calculations.— 

Applying the simple criterion described above, Ca and Sr appear to be barely stable in 
their surface location with respect to the bulk one,^^ as reflected in the A values collected 
in Table [H which are close to Aq for these doped ^He systems. This borderline character for 
the solvation properties of these impurities implies that detailed calculations are required to 
help to understand the results of spectroscopic studies on alkalineearth- doped He droplets. 
In particular, high-quality impurity- He pair interaction potentials are required since even 
relatively small inaccuracies in these potentials, which are often not known with a sufficient 
precision, may yield wrong results. 

We present here a systematic study for helium drops made of each isotope, having a 
number of atoms large enough to make them useful for the discussion of experiments on 
laser-induced fluorescence (LIF) or beam depletion (BD) spectroscopy or for the discus- 
sion of other physical phenomena involving these systems, such as interatomic Coulombic 
decay.— We also discuss the dependence of the structural properties on the cluster size. 
Some of this information is also experimentally available.— After a brief explanation of the 
experiment, results are presented for strontium on helium nanodroplets that further support 
the calculations. 

This work is organized as follows. In Sec. [Tll we briefly describe the density functional 
plus alkaline earth-He potential approach employed here, as well as some technical details. 
Doped drops calculations are presented and discussed in Sec. Illlt while the experimental 
results are discussed in Sec. HVt and an outlook is presented in Sec. |Vl 
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II. DENSITY-FUNCTIONAL DESCRIPTION OF HELIUM NANODROPLETS 



Since the pioneering work of Stringari and coworkers,— density functional (DF) theory 
has been used in many studies on hquid hehum in confined geometries and found to provide 
a quite accurate description of the properties of inhomogeneous hquid He (see e.g. Ref. [j 
and references therein). 

The starting point is to write the energy of the system as a functional of the He particle 
density p: 

E[p] = j drS{p) + J dr'p{r')VAke~He{\r-r'\) , (1) 

where S{p) is the He energy density per unit volume, and VAke-He is the alkaline earth- helium 
pair potential. The impurity is thus treated as a fixed external potential. Addressing the 
lightest alkaline earth. Be, for which a fairly recent Be-He is available," would likely require 
to treat this atom as a quantum particle instead of as an external potential.— 

For ^He we have used the Orsay-Trento functional,— and for ^He the one described in 



Ref. |28| and references therein. These functionals have been used in our previous work on 
helium drops doped with alkali atoms^*^ as well as in many other theoretical works. The 
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results discussed in the following have been mostly obtained using the potentials of Ref. 
(Ca, Sr and Ba), and of Ref. 30 (Mg, for which the pair potentials of Refs. 26 andjsol are 
similar). For Ca, we have also tested other potentials available in the literature,— >^>^ as 
well as the unpublished potential of Meyer— we had employed in our previous work.— 
Fig. [1] shows the pair potentials used in this work. From this figure, one may anticipate 



that Ca@^Heiv drops described using the potentia 
the same drops described with the potential of Ref. 
potentials of Refs. 



of Ref. 



26l and 



321 display deeper dimples than 
29| . We want to point out that the Ca-He 
and should yield equivalent 
is more attractive, causing 
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3d are very similar to that of Ref. 
results. Contrarily, we have found that the potential of Ref. 
the Ca atom to be drawn to the center of the ^Hctv drop, in contrast with the experimental 
findings.— 

For a number of helium atoms in the drop, we have solved the Euler-Lagrange equation 
which results from the variation of E[p] at constant A^: 



T \- VAke-He 

5p 



1^, (2) 
where p is the helium chemical potential, whose value is determined self-consistently by 
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imposing the auxiliary condition / drp{r) = N during the iterative minimization. 

When the impurity resides off center (as in the case of a dimple structure), the system is 
axially symmetric. Despite this symmetry, we have solved Eq. ([2]) 

in Cartesian coordinates because this allows us to use fast Fourier transform techniques'^ 
to efficiently compute the convolution integrals entering the definition of S{p), that is, the 
mean field helium potential and the coarse-grained density needed to evaluate the correlation 
term in the density functional.— We have found this procedure to be faster and more accurate 
than convoluting by direct integration using cylindrical coordinates. 

We have used an imaginary time method'^'^ to solve Eq. ([2]), after having discretized it 
using 13-point formulas for the spatial derivatives. The mesh used to discretize p in space 
is chosen so that the results are stable against small changes of the mesh step. 



III. THEORETICAL RESULTS 



We start a typical calculation by placing the impurity close to the surface of the He 
droplet. Depending on the studied impurity and/or the He isotope, during the functional 
minimization, the alkaline earth atom is either driven to the interior of the droplet,— or it 
remains trapped in a more or less pronounced dimple on its surface. 

In the case of ^He, we find that for all of the alkaline earth atoms investigated, the stable 
state is always the one where the impurity is in the center of the cluster. This is consistent 
with the associated large A values, see Table [H Figure [2] shows the density profiles for 
MgOSHejv, Ca@3He^, SrO^He^, and Ba@3HeAr for N = 300, 500, 1000, 2000, 3000, and 
5000. For CaQ^Hesooo, we also show the profile obtained with the pair potential of Ref. |32 



(dotted line). Several solvation shells are clearly visible. The number of ^He atoms below 
the first solvation peak for the = 5000 drop is about 19 for Mg, 22 for Ca, 26 for Sr, 
and 27 for Ba. The differences in the location and height of the first solvation peak are a 
simple consequence of the different depth and equilibrium distance of the corresponding pair 
potentials. It is interesting to see the building up of the drop structure around the impurity 
that, as it is wellknown,^ only causes a large but localized effect on the drop structure. 

The bottom panel of Fig. [3] shows the corresponding solvation energies, defined as the 
energy differences 

S^iAke) = E{Ake@^HeN) - E^^Hcn) , (3) 
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with an equivalent definition for ^He drops. Tlie more attractive Ca-He pair potential of 



Ref. 



32 



yields, on average, a solvation energy about 13 K larger as compared with that 



obtained with the pair potential of Ref. |29|, despite the fact that the density profiles look 
fairly similar; see Fig. [2J 

In the case of Ca and Sr atoms in ^He drops, whose A values are close to the threshold 
for solvation Aq (see Table H]), we have found that, for both dopants, the minimum energy 
configuration is a dimple state at the surface, although the energy difference between the 
surface and the solvated states is fairly small for both dopants. For Ca@'*He3oo, this difference 
is 3.4 K using Meyer's potential,— and 12.0 K using that of Ref. |29|. The homologous result 
for Sr@'*He3oo is 22.7 K. These energy differences have to be compared with the total energy 
of the ^Heaoo drop, which is about —1384 K. 

We have also confirmed by DF calculations the surface state of Ba@'*HeAr and the sol- 
vated state of Mg@''He7v, both suggested by the corresponding A values in Table [B This is 
illustrated in Fig. H] for Mg, and in Fig. [5] for Ca, Sr and Ba. The dimple depth ^, defined 
as the difference between the position of the dividing surface at p = pfe/2 -where ph is the 
bulk liquid density- with and without impurity, respectively, is shown in Fig. [6] as a func- 
tion of A^. The structure of the dimple is different for different alkaline earth atoms, being 
shallower for Ba and more pronounced for Ca. We recall that the ^ values for Na@^He2ooo 
and Na@^He2ooo are 4.5 and 2.1 A, respectively.- The dimple depths for alkaline are thus 
much smaller than for alkaline earth atoms, as also indicated by LIF experiments.-'^'^'^ 
The dependence of the the dimple depth with the alkaline earth atom size, characterized 
by the radial expectation value RAke of the valence electrons,— is shown in Fig. [71 This 
figure is consistent with the increasing bulk-to-surface ratio of the line shifts as the size of 
the dopant atom increases.— 

The 'solvation' energies for these alkaline earth atoms in ^He drops are displayed in the 
top panel of Fig. [31 As in the case of ^He drops discussed before, the stronger the Ake- 
He pair potential (see Fig. [T]), the more negative SN^Ake). In the case of Ca@^HeAr and 
Sr@^HeAr, the energies are very similar, and so are the dimple depths shown in Fig. [6l It is 
worth seeing the different behavior of Sn as a function of A^ for each helium isotope. In the 
case of ^He, once the first 2-3 solvation shells are fully developed, S'a? quickly saturates, and 
for this reason it changes only by 12 % (Ca) and 17 % (Sr) from A^ = 300 to A^ = 5000. For 
the same reason, spectroscopic shifts are expected to be A^ independent for drops made of 
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more than a few hundred ^He atoms. When the impurity is at the surface, sizeable curvature 
effects appear even for a few thousand atoms drops. This shows up not only in the change 
of Sm, which is about 22 % for Ca, and 24 % for Sr in the same N range as before, but 
also in the spectroscopic shifts, that still depend on N below ~ 3000 (see e.g. Ref. llTI ). 
This illustrates the need of large drops for carrying out spectroscopic shift calculations to 
attempt a detailed comparison with experiments. 



IV. EXPERIMENTAL RESULTS 

To support the DF calculations, the 5s5p ^P^ ^ ^Sq transition of strontium on nan- 
odroplets made of either helium isotope has been experimentally investigated. Although 
calcium appears to be most favorable, we are so far restricted to excitation spectra of stron- 
tium attached to helium droplets because of the limited tuning range of our lasers. Calcium 
will be addressed in a future experiment. The experiments where performed in a helium 
droplet machine applying laser-induced fluorescence, as well as beam depletion and photo 
ionization (PI) spectroscopy. A detailed description of the experimental setup is presented 
elsewhere.— Modifications include a new droplet source to reach the lower temperatures 
needed for generating ^He droplets.-^ In short, gas of either helium isotope is expanded un- 
der supersonic conditions from a nozzle, forming a beam of droplets traveling freely in high 
vacuum. The helium stagnation pressure in the droplet source is 20 bar, and a nozzle of 
5 /im diameter has been used. The nozzle temperature has been stabilized to 12 K and 15 K 
to form ^He and ^He droplets, respectively. These conditions result in an average droplet 
size of ~ 5000 helium atoms.- 

The droplets are doped downstream using the pick-up technique: in a heated scattering 
cell, an appropriate vapor pressure of strontium is established so that droplets pick up 
one single atom on average when passing the cell. LIF as well as PI and BD absorption 
spectra of the doped droplet beam can be recorded upon electronic excitation using a pulsed 
nanosecond dye laser. LIF is recorded with a photo multiplier tube. In the case of PI, the 
photons of an excimer laser ionize the excited atoms in a one photon step. The ions are 
afterwards detected by a channeltron. For the beam depletion measurement, a Langmuir- 
Taylor surface ionization detector has been used.— 

We show here only the results obtained using LIF because a much better signal-to-noise 
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ratio was achieved when compared to the PI and BD spectra for strontium doped clusters. 
In the case of PI, the reason probably is the tendency of the just formed strontium ions 
not to desorb from the droplet like e.g. alkali atoms do. Since the detection efficiency of 
our PI detector is considerably decreased for high masses, the detection of the ion+droplet 
complex is small. A decreased desorption mechanism also diminishes the sensitivity of BD 
techniques. However, the PI/BD measurements give identical results when compared to the 
LIE spectra. 

Figure [8] shows the measured spectra of the 5s5p ^P° ^ ^Sq transition of strontium 
atoms on droplets of '^He/^He compared to that in bulk ^He.^^ All three spectra show a broad 
asymmetric line, blue shifted from the atomic gas-phase absorption. The differences of the 
shifts for ^He drops and bulk ^He immediately confirms the surface location of the strontium 
atoms. In the case of bulk ^He, the absorption is far more blue shifted and the width is 
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considerably wider. The shift can be explained within the bubble model, see e.g. Refs. |1 
and references therein, and results from repulsion of the helium environment against spatial 
enlargement of the electronic distribution of the excited state. The shift in bulk helium is 
larger than in droplets because the dopant is completely surrounded by helium, whereas it 
is not when it is located at the surface of drops. 

Table HTl summarizes the experimentally determined shifts of the first electronic transition 
of strontium and calcium in helium droplets, as well as the measurements in bulk helium 
for both isotopes. As compared to ''He drops, the absorption maximum in the case of ^He 
drops is shifted 60 cm~^ further to the blue, and the width increases from 180 to 220 cm^^. 

At first glance, it is not obvious from the recorded spectra in ^He drops whether the 
strontium atom is in a surface state, or it is solvated inside the droplets. It is worth men- 
tioning that Morowaki et al. performed similar measurements in bulk helium.— They have 
compared the absorption spectra of the 4s4p ^P° ^ ^Sq transition of calcium in bulk 
^He and '^He, and have found a much smaller blue shift in the case of ^He (about 55%, see 
Table Hill . which could again be explained within the bubble model -the reduced shift just 
results from the lower density of liquid ^He. A similar quantitative effect should be expected 
for strontium, especially in view of the reported DF calculations. 

Consistent with this expectation is that in our experiments, the measured shift of Sr 
in ^He droplets, 140 cm~^, is about a 58% of the value corresponding to Sr in bulk ^He, 
240 cm~^.-^ We can safely argue that the shift determined in ^He droplets should sensibly 
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coincide with the expected value for bulk ^He, indicating complete solvation of strontium 
atoms in ^He droplets, as predicted by DF calculations. 

V. SUMMARY AND OUTLOOK 

In this work, we have presented detailed results for the structure and energetics of helium 
drops doped with Mg, Ca, Sr, and Ba alkaline earth atoms. We have found that these atoms 
are solvated in the case of ^He drops and reside in surface dimples in the case of ^He drops, 
with the sole exception of Mg@^He7v, which is also solvated. This yields a fairly complete 
physical picture, from the theoretical viewpoint, of the structure and energetics of helium 
drops doped with alkaline earth atoms. The experimental spectrum of strontium atoms in 
^He and ^He droplets confirms the DF calculations. Moreover, since the spectroscopic shift is 
sensitive to the shape/depth of the surface dimple, a comparison between experimental and 
calculated line shifts could provide a sensible test on the accuracy of available pair potentials. 
We want to stress again that accurate pair potentials are needed to quantitatively reproduce 
the experimental results, especially when the solvation properties of the impurity are such 
that they yield values of A close to the threshold value Aq. 

The different solvation behavior of the heavier alkaline earth atoms in ^He and ^He drops, 
offers the unique possibility of using them to study mixed drops at very low temperatures, in 
particular the ^He-'^He interface. It is known that below the tricritical point at ~ 0.87 K,— 
^He has a limited solubility in ^He, segregating for concentrations larger than a critical 
value. This segregation also appears in mixed droplets,— i^i^ producing a shell structure in 
which a core, essentially made of ^He atoms, is coated by ^He that is hardly dissolved into 
the ^He core, even when the number of ^He atoms is very large. Due to this particular 
structure, that pertains to medium to large size droplets, strongly attractive impurities reside 
in the *^He core, being very little affected by the outer ^He shell, whereas weakly attractive 
impurities, like alkali atoms, should still reside in the surface of the droplet, irrespective 
of the existence of the ^He core. Contrarily, Ca, Sr and Ba impurities would be sunk into 
the fermionic component up to reaching the ^He-^He interface if the appropriate number 
of atoms of each isotope is chosen. This will offer the possibility of studying the ^He-^He 
interface, and a richer alkaline earth atom environment. We are at present generalizing the 
DF approach we have used in the past^^ to address this more demanding and promising new 
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aspect of the physics of doped hehum droplets. On the experimental side, calcium spectra 
will be accessible in forthcoming experiments. We want to point out that mixed droplets 
doped with alkali atoms have been already detected in our previous experiments,— and that 
systematic experiments on alkaline earth doped mixed droplets will be performed in the 
future. 
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TABLE I: A parameter for the alkaline earth atoms and pair potentials used in this work. 
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TABLE II: Experimental shifts of the first electronic transition of Ca and Sr atoms in bulk helium 
as well as in drops. The values for Sr@HeAr are from this work. Previous experiments, carried out 
only for Sr@^HeAr, showed the same shifts 
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FIG. 1: Alkaline earth-He pair 
of doped helium drops: (1) Ref. 



potentials used in this work to obtain the ground state structure 



30; (2) Ref. 
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FIG. 2: Density profiles for ^HeTv drops doped with Mg, Ca, Sr, and Ba, for N = 300, 500, 1000, 
2000, 3000, and 5000. The dotted line in the Ca panel corresponds to CaQ^Hesooo calculated with 



the pair potential of Ref. 



potentials of Ref. 



32l . Drops doped with Ca, Sr and Ba have been calculated using the pair 



29! , and drops doped with Mg, using the pair potential of Ref. 
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FIG. 3: Top panel: solvation energies (K) for doped 



30; (2) from Ref. 



^HeAT drops. Resu 



2g|; (3) from Ref. 



following pair potentials: (1) from Ref. 
as top panel for doped '^HeTv drops. The lines are drawn to guide the eye 



ts obtained using the 



32l . Bottom panel: same 
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FIG. 5: Equidensity lines on a symmetry plane for ^Hejv drops with N =300 (left panels) and 
1000 (right panels) doped with Ca, Sr and Ba. The lines span the surface region between 0.9pb 
and O.lph in O.lph steps, where pb is the bulk liquid density 0.0218 A^^. The cross indicates the 
location of the alkaline earth atom in the dimple. Results obtained using the pair potentials of 
Ref. 29. 
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FIG. 6: Depth of the dimples (^) created m ^HeA? drops obtained using the fohowing pair potentials: 



(2) from Ref. 



23 for Ba (diamonds), Sr (circles), and Ca (solid dots) atoms; (3) from Ref. 



32 for 



Ca (squares). The lines are drawn to guide the eye. 
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FIG. 7: Depth of the dimples (^) created in ^Hesooo drops 



of the atomic size R^f^^, using the pair potentials of Ref. 



29 



ay Ba, Sr and Ca atoms, as a function 
The line is drawn to guide the eye. 
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